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BACKGROUND:  Hypothermia,  which  can  result  from  tissue  hypoperfusion,  body  exposure,  and 
transfusion  of  cold  resuscitation  fluids,  is  a  major  factor  contributing  to  coagulopathy  of  trauma 
and  surgery.  Despite  considerable  efforts,  the  mechanisms  of  hypothermia-induced  blood  coag¬ 
ulation  impairment  have  not  been  fully  understood.  We  introduce  a  kinetic  modeling  approach  to 
investigate  the  effects  of  hypothermia  on  thrombin  generation. 

METHODS:  We  extended  a  validated  computational  model  to  predict  and  analyze  the  impact 
of  low  temperatures  (with  or  without  concomitant  blood  dilution)  on  thrombin  generation  and 
its  quantitative  parameters.  The  computational  model  reflects  the  existing  knowledge  about 
the  mechanistic  details  of  thrombin  generation  biochemistry.  We  performed  the  analysis  for  an 
“average”  subject,  as  well  as  for  472  subjects  in  the  control  group  of  the  Leiden  Thrombophilia 
Study. 

RESULTS:  We  computed  and  analyzed  thousands  of  kinetic  curves  characterizing  the  generation 
of  thrombin  and  the  formation  of  the  thrombin-antithrombin  complex  (TAT).  In  all  simulations, 
hypothermia  in  the  temperature  interval  31°C  to  36°C  progressively  slowed  down  thrombin  gen¬ 
eration,  as  reflected  by  clotting  time,  thrombin  peak  time,  and  prothrombin  time,  which  increased 
in  all  subjects  (P  <  10-5).  Maximum  slope  of  the  thrombin  curve  was  progressively  decreased, 
and  the  area  under  the  thrombin  curve  was  increased  in  hypothermia  (P  <  10-5);  thrombin  peak 
height  remained  practically  unaffected.  TAT  formation  was  noticeably  delayed  (P  <  10-5),  but  the 
final  TAT  levels  were  not  significantly  affected.  Hypothermia-induced  fold  changes  in  the  affected 
thrombin  generation  parameters  were  larger  for  lower  temperatures,  but  were  practically  inde¬ 
pendent  of  the  parameter  itself  and  of  the  subjects’  clotting  factor  composition,  despite  sub¬ 
stantial  variability  in  the  subject  group.  Hypothermia  and  blood  dilution  acted  additively  on  the 
thrombin  generation  parameters. 

CONCLUSIONS:  We  developed  a  general  computational  strategy  that  can  be  used  to  simulate 
the  effects  of  changing  temperature  on  the  kinetics  of  biochemical  systems  and  applied  this 
strategy  to  analyze  the  effects  of  hypothermia  on  thrombin  generation.  We  found  that  thrombin 
generation  can  be  noticeably  impaired  in  subjects  with  different  blood  plasma  composition  even 
in  moderate  hypothermia.  Our  work  provides  mechanistic  support  to  the  notion  that  thrombin 
generation  impairment  may  be  a  key  factor  in  coagulopathy  induced  by  hypothermia  and  compli¬ 
cated  by  blood  plasma  dilution.  (Anesth  Analg  2013;117:565-74) 
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Hypothermia,  a  state  of  abnormally  low  body  core 
temperature  (<36°C  in  trauma  patients),  is  one  of 
the  major  factors  contributing  to  coagulopathy 
of  trauma.1-3  Heavy  bleeding  may  result  in  a  decreased 
blood  volume  and  associated  tissue  hypoperfusion,  which, 
in  turn,  leads  to  insufficient  heat  generation  in  the  tissues 
due  to  a  reduced  rate  of  oxygen-dependent  adenosine  tri¬ 
phosphate  synthesis.4  Transfusion  of  blood  products  and 
resuscitation  fluids,  which  are  stored  at  low  temperatures 
or  even  in  a  frozen  state,  further  contribute  to  the  develop¬ 
ment  of  hypothermia  and  coagulopathy.2-3  Body  rewarming 
techniques  are  generally  expected  to  mitigate  the  adverse 
effects  of  hypothermia,1-3  but  they  might  require  time  to  take 
effect  when  rapid  action  is  needed  for  a  patient's  survival. 
Moreover,  in  trauma  situations,  hypothermia-induced 
coagulopathy  may  be  complicated  by  other  factors  (such  as 
blood  dilution),5  which  cannot  be  corrected  by  rewarming. 
To  develop  new,  broad-spectrum  therapeutic  strategies  for 
reversing  traumatic  and  surgical  coagulopathy,6  it  is  thus 
necessary  to  investigate  the  effects  of  hypothermia  on  blood 
clotting  at  the  molecular  level. 

The  impairment  of  blood  clotting  under  hypothermic 
conditions  is  well  recognized.5  However,  the  mechanisms  of 
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Reduced  Temperature  Effects  on  Thrombin  Generation 


such  impairment  have  not  been  firmly  established.  Some  of 
the  existing  experimental  evidence  indicates  that  hypother¬ 
mia  noticeably  impacts  the  generation  of  thrombin,  a  central 
enzymatic  component  of  the  biochemical  blood  coagulation 
network.7-10  Yet,  other  reports  suggest  that  thrombin  gen¬ 
eration  biochemistry  is  practically  unaffected  by  mild  and 
moderate  hypothermia  and  emphasize  the  role  of  hypo¬ 
thermia-induced  effects  that  are  not  attributed  to  thrombin 
generation  deficiency.11-14  The  use  of  different  experimental 
systems  and  conditions  by  different  laboratories,  which  can 
quantitatively  impact  thrombin  generation  kinetics,15,16  has 
hampered  definitive  conclusions. 

The  ambiguity  regarding  the  effects  of  temperature  on 
thrombin  generation  suggests  that  experimental  results 
should  be  analyzed  and  interpreted  using  new,  complemen¬ 
tary  approaches.15  Here,  we  apply  computational  methods 
to  help  understand  and  predict  temperature  modulation  of 
thrombin  kinetic  curves.  This  approach,  based  on  extending 
the  validated  Hockin-Mann  computational  model  of  throm¬ 
bin  generation,17,18  was  motivated  by  recent  applications  of 
this  model  to  study  the  impairment  of  thrombin  generation 
in  dilution-induced  coagulopathy19  and  to  investigate  the 
therapeutic  effects  of  factor  Vila15  and  prothrombin  complex 
concentrates20  on  thrombin  accumulation.  Moreover,  the 
Hockin-Mann  model  has  been  of  considerable  utility  in  the 
analysis  of  other  clinically  relevant  blood  clotting  problems, 
such  as  causes  for  thromboembolic  complications21,22  and 
the  action  of  anticoagulants.23,24  The  model  reflects  the  cur¬ 
rent  mechanistic  view,  albeit  simplified,  of  the  biochemical 
thrombin  generation  network  in  human  blood  plasma,  and 
captures  the  typical  kinetics  of  thrombin  generation  in  vitro. 

Here,  we  extend  the  capabilities  of  the  Hockin-Mann 
model  and  introduce  a  rigorous  and  general  computational 
strategy  to  model  the  effects  of  changing  temperature  on 
the  kinetics  of  biomolecular  systems.  We  use  this  strategy  to 
analyze  the  temperature  dependence  of  thrombin  genera¬ 
tion  kinetics.  The  contributions  of  this  work  are  2-fold.  First, 
we  characterize  the  effects  of  hypothermia  (alone  and  in 
combination  with  blood  dilution)  on  the  5  standard  quanti¬ 
tative  thrombin  generation  parameters,  which  has  not  been 
previously  done.  Second,  we  establish  how  the  effects  of 
hypothermia  on  thrombin  generation  are  impacted  by  inter¬ 
subject  variability. 

METHODS 
Study  Group 

Because  this  work  is  a  computational  study  that  uses  pre¬ 
viously  obtained  experimental  data,  no  IRB  approval  was 
necessary.  We  performed  computational  analyses  for  a 
subject  characterized  by  average  values  of  the  concentra¬ 
tions  of  blood  coagulation  components  in  normal  human 
blood  plasma  (the  "average"  subject),15,18  as  well  as  for 
472  individual  subjects  in  the  control  group  of  the  Leiden 
Thrombophilia  Study  (LETS).25 

Computational  Modeling 

The  Hockin-Mann  model  is  a  system  of  mathematical  equa¬ 
tions  where  the  unknown  variables  are  the  concentrations 
of  the  biochemical  species  in  the  thrombin  generation  net¬ 
work.17,18  The  equations  reflect  the  chemical  reactions  in  the 


network  (listed  in  Table  Al,  see  Appendix)  and  link  the  con¬ 
centrations  with  their  rates  of  change.  The  solutions  of  the 
equations  represent  the  temporal  dynamics  of  the  concen¬ 
tration  variables.  The  model  is  parameterized  by  the  values 
of  44  constants  (termed  kinetic  constants)  that  define  the 
rates  of  the  reactions  in  the  thrombin  generation  network. 
The  input  of  the  model  is  the  initial  concentrations  of  the 
proteins  that  participate  in  the  reactions,  and  the  model's 
output  is  the  resulting  kinetic  trajectories  for  the  protein 
concentrations  computed  over  a  specified  time  interval  (50 
minutes).  The  default  model  input  is  the  average  coagula¬ 
tion  factor  concentrations  in  normal  human  plasma,  and 
this  input  can  be  changed  to  represent  alternative  blood 
compositions,  such  as  those  in  the  LETS  subjects.  The  LETS 
dataset  contains  relative  levels  of  blood  coagulation  factors 
for  each  of  the  subjects.  We  rescaled  them  using  the  known 
average  concentration  values  in  absolute  units  and  used  the 
resulting  concentrations  as  the  model  inputs  to  simulate 
thrombin  generation  in  individual  subjects.  In  our  simu¬ 
lations  for  diluted  blood  plasma,  the  default  inputs  were 
modified  according  to  the  dilution  scenarios.19,20 

The  extension  of  the  Hockin-Mann  model  proposed  in 
this  work  concerns  the  effects  of  changing  temperature  on 
biochemical  kinetics.  Because  accurate  values  of  the  param¬ 
eters  (termed  temperature  coefficients  [TCs])  that  gov¬ 
ern  those  effects  have  not  been  determined,  our  approach 
is  based  on  random  sampling  of  TC  sets  and  using  the 
Hockin-Mann  model  to  compute  a  thrombin  curve  (i.e.,  a 
kinetic  trajectory  for  the  concentration  of  thrombin)  for  each 
of  such  sets.  Groups  of  thrombin  curves  computed  for  many 
sampled  TC  sets  define  the  predicted  kinetic  curve  ranges 
(or,  "corridors")  to  which  the  true  hypothermic  thrombin 
curves  are  expected  to  belong. 

A  typical  procedure  for  modeling  temperature  depen¬ 
dence  can  be  described  as  follows.  For  any  given  hypo¬ 
thermic  temperature  T2  (all  temperatures  are  expressed 
in  degrees  Celsius)  and  each  of  the  44  kinetic  constants 
in  the  model,  we  generated  a  number,  M,  of  random  TCs 
(see  Appendix  for  details  on  TC  sampling).26,27  (The  use  of 
TCs  is  grounded  in  physical  chemistry,  being  an  approxi¬ 
mation  of  the  well-known  Arrhenius  equation.26,28)  M  was 
equal  to  5000  for  most  simulations.  As  a  result  of  such  ran¬ 
dom  sampling,  we  obtained  M  sets,  each  set  containing  44 
TC  values.  In  every  set,  each  TC  was  used  to  calculate  the 
hypothermic  value  of  1  of  the  44  kinetic  constants  as  fol- 
lows:  k2  =  /q  x  y  2  1  ,  where  k2  and  k2  denote  the  values  of 

the  kinetic  constant  at  temperatures  Tj  and  T2,  respectively; 
y  denotes  the  TC;  and  Tj  =  37°C  represents  the  normal  body 
temperature.26,27  For  each  of  the  M  sets,  the  equations  con¬ 
stituting  the  Hockin-Mann  model  were  solved  numerically, 
and  a  thrombin  trajectory  was  generated.  The  M  resulting 
trajectories  were  analyzed  using  statistical  methods,  and 
predictions  of  thrombin  generation  at  the  hypothermic  tem¬ 
perature  T2  were  obtained  from  those  analyses.  Additional 
details  about  our  computational  procedures  are  given  in  the 
Appendix. 

Statistical  Analysis 

Groups  of  simulated  thrombin  curves  and  thrombin- 
antithrombin  complex  (TAT)  curves  were  analyzed  using 
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descriptive  statistics.  (These  curve  groups  were  generated 
by  TC  randomization,  as  described  in  the  previous  subsec¬ 
tion.)  The  thrombin  or  TAT  raw-output  kinetic  curves  from 
the  computational  model  were  discretized  into  500  evenly 
spaced  points.  Each  of  these  points  corresponded  to  a  given 
time  moment  within  the  default  simulation  time  interval. 
For  each  of  these  time  moments,  we  analyzed  the  statistical 
sample  of  thrombin  or  TAT  levels  constituted  by  the  corre¬ 
sponding  curve  values  from  every  curve  in  the  group.  For 
each  of  the  resulting  500  samples,  we  calculated  the  range 
(i.e.,  maximum  variation),  interquartile  range,  and  the 
median,  which  were  plotted  against  time  to  represent  their 
temporal  changes  (Figs.  1  and  2). 

For  groups  of  model-generated  thrombin  and  TAT 
curves,  we  calculated  and  statistically  analyzed  quantita¬ 
tive  parameters,  such  as  clotting  time  (CT)  and  thrombin 
peak  height  (PH)  (see  Results).  In  these  analyses,  the  curves 
in  a  group  either  were  generated  for  a  given  subject  via  TC 
randomization  or  were  generated  for  every  subject  in  the 
LETS  subject  group  in  the  absence  of  TC  randomization.  For 


each  curve  in  the  group,  the  quantitative  parameters  were 
calculated.  The  parameter  groups  constituted  the  statistical 
samples  for  which  we  calculated  interquartile  ranges  and 
means  and  performed  sample  comparisons  using  nonpara- 
metric  hypothesis  testing.  Kinetic  curve  parameter  samples 
were  compared  using  Wilcoxon  rank  sum  test  when  TCs 
were  randomized  (in  the  computations  for  the  "average" 
subject)  or  Wilcoxon  signed  rank  test  when  no  TC  random¬ 
ization  was  performed  (in  the  computations  for  the  LETS 
subjects).  The  resulting  P  values  were  Bonferroni-corrected 
to  account  for  multiple  comparisons  of  a  thrombin  param¬ 
eter  distribution  for  the  normal  temperature  with  the  dis¬ 
tributions  of  the  same  parameter  for  different  hypothermic 
temperatures  (for  a  fixed  dilution  scenario). 

RESULTS 

Temperature  Modulation  of  Thrombin  Curves 

For  each  of  the  considered  hypothermic  temperatures 
(31  °C,  32°C,  33°C,  34°C,  35°C,  and  36°C),  we  independently 


Figure  1.  Thrombin  generation  under 
normal  and  hypothermic  conditions. 
TAT  designates  thrombin-antithrombin 
complexes,  and  brackets  designate  con¬ 
centration.  A  and  B,  Model-generated 
thrombin  curves  for  the  “average”  sub¬ 
ject  under  normal  temperature  (blue) 
and  hypothermic  temperature  (red) 
conditions.  Each  red  curve  corresponds 
to  1  of  30  randomly  generated  sets  of 
temperature  coefficients  (TCs).  C  and 
D,  Model-generated  thrombin  curves  for 
the  average  subject  under  normal  tem¬ 
perature  (blue)  and  median  thrombin 
curves  for  simulated  hypothermia.  (Five 
thousand  random  TC  sets  were  gener¬ 
ated  independently  for  each  hypothermic 
temperature.  The  resulting  thrombin  and 
TAT  curves  were  used  to  calculate  the 
median  thrombin  and  TAT  values  at  each 
of  the  time  points  that  discretize  the 
considered  time  interval.)  E,  Normalized 
data  on  hypothermia  induced  in  vivo  in  a 
porcine  model;  data  taken  from  Martini 
et  a  I.29 
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Figure  2.  Simulated  thrombin  generation  ranges  for  hypothermic  temperatures  in  the  “average"  subject.  A,  Thrombin  curves;  B,  Thrombin- 
antithrombin  complex  (TAT)  curves.  Blue  curves  correspond  to  normal  temperature.  Solid  black,  green,  and  red  curves  correspond  to  the 
range,  interquartile  range,  and  median,  respectively,  calculated  from  the  hypothermic  curves  generated  for  5000  random  temperature  coef¬ 
ficient  (TC)  sets.  The  dashed  black  lines  correspond  to  hypothermic  curves  calculated  using  a  model  for  which  all  TCs  were  assigned  the  value 
2.5  (no  randomization).  All  hypothermic  thrombin  and  TAT  curves  were  generated  for  33°C. 


generated  5000  random  sets  of  TCs,  and  for  each  TC  set  we 
calculated  thrombin  and  TAT  curves  for  the  average  subject. 
Because  the  generated  curves  depended  on  the  TCs,  there 
were  variations  between  the  curves  within  each  curve  group 
(with  1  group  corresponding  to  1  hypothermic  tempera¬ 
ture).  Yet,  for  all  the  model-generated  hypothermic  thrombin 
curves,  CT  (time  to  10  nM  thrombin)  and  thrombin  peak  time 
(PT)  were  larger  than  the  average  subject  values  at  normal 
temperature,  i.e.,  37°C.  The  same  result  held  for  the  50%  acti¬ 
vation  time  (the  time  it  takes  the  curve  to  reach  50%  of  its  final 
plateau  value,  a  standard  measure  of  timing  for  the  activa¬ 
tion  of  biological  processes)30  for  all  hypothermic  TAT  curves 
(Fig.  1,  A-D).  These  results  are  consistent  with  a  number  of 
works  reporting  hypothermia-delayed  clotting  measured 
by  thromboelastography  and  related  techniques  in  human 
blood,7-9-31'32  as  well  as  with  TAT  measurements  for  a  porcine 
model  of  hypothermia-induced  coagulopathy  (Fig.  IE). 

Localization  of  Predicted  Thrombin  Curves  Under 
Hypothermic  Conditions 

We  investigated  how  well  the  median  curves  describe  the 
behavior  of  groups  of  model-generated  thrombin  and  TAT 
curves  at  each  hypothermic  temperature.  For  every  group 
of  5000  curves  described  in  the  previous  subsection,  we 
calculated  the  curves  corresponding  to  the  range  and  inter¬ 
quartile  range.  The  plots  for  33°C  are  shown  in  Figure  2; 
the  plots  for  all  other  considered  hypothermic  tempera¬ 
tures  were  similar  (results  not  shown).  While  the  calculated 
ranges  for  the  hypothermic  curves  could  have  substantial 
width  (due  to  the  sensitivity  to  outliers),  the  interquartile 
ranges  for  all  predicted  hypothermic  thrombin  and  TAT 
curves  were  quite  narrow  (Fig.  2).  Therefore,  if  the  sample 
medians  were  estimated  accurately,  they  could  serve  as  a 
proxy  for  the  behavior  of  the  hypothermic  curves  for  the 
"true"  (unknown)  TCs. 

To  assess  the  accuracy  of  estimation  of  the  median 
curves,  we  repeated  all  described  simulations  for  a  larger 
sample  size  (10,000  random  sets  of  TCs,  generated  indepen¬ 
dently  for  every  considered  hypothermic  temperature).  The 
resulting  curves  for  the  medians  and  interquartile  ranges 
practically  coincided,  and  the  curves  for  the  ranges  nearly 
coincided,  with  the  corresponding  curves  for  sample  size 


5000  (results  not  shown).  Thus,  sample  size  5000  was  used 
in  all  subsequent  analyses  involving  TC  randomization. 

Because  generation  and  analysis  of  thousands  of  throm¬ 
bin  curves  are  computationally  expensive,  we  attempted  to 
simplify  the  method.  We  hypothesized  that  hypothermic 
thrombin  and  TAT  kinetic  curves  calculated  for  the  aver¬ 
age  subject  with  each  TC  equal  to  2.5  (middle  of  the  sam¬ 
pling  interval;  see  Appendix)  could  serve  as  a  good  estimate 
for  the  median  curves  calculated  with  TC  randomization. 
Comparison  of  the  curves  showed  a  very  good  correspon¬ 
dence  for  all  considered  temperatures  (the  case  of  33°C  is 
shown  in  Fig.  2). 

Combined  Effects  of  Temperature  and  Dilution 
on  Standard  Quantitative  Parameters  of 
Thrombin  Generation 

For  the  generated  groups  of  5000  thrombin  curves  for  each 
hypothermic  temperature,  we  calculated  the  5  standard 
thrombin  curve  parameters15-19:  CT,  PT,  maximum  slope 
of  the  thrombin  curve  (MS),  PH,  and  the  area  under  the 
thrombin  curve  (AUC).  Moreover,  we  used  the  same  strat¬ 
egy  to  generate  thrombin  curves  and  calculate  the  quanti¬ 
tative  parameters  for  each  hypothermic  temperature  when 
the  blood  plasma  in  the  average  subject  was  diluted  3-  or 
5-fold.  The  resulting  parameter  distributions  are  shown  in 
Figure  3,  A  to  E.  For  each  considered  temperature,  the  medi¬ 
ans  and  interquartile  ranges  of  the  parameter  distributions 
followed  a  pattern  that  was  consistent  with  our  previous 
results  for  normal  temperature19:  dilution  increased  CT  and 
PT,  decreased  MS  and  PH,  and  left  AUC  almost  unchanged; 
these  effects  were  more  pronounced  for  5-fold  than  for 
3-fold  dilution.  As  the  temperature  decreased,  CT  and  PT 
monotonically  increased,  MS  monotonically  decreased,  and 
PH  remained  almost  unchanged  (Fig.  3,  A-D).  Remarkably, 
AUC  demonstrated  a  robust  increase  for  decreasing  tem¬ 
perature  (Fig.  3E). 

To  model  the  effects  of  hypothermia  on  prothrombin 
times  (typically  reported  in  experimental  studies),  we  per¬ 
formed  similar  calculations  for  CT  when  thrombin  genera¬ 
tion  was  activated  by  tissue  factor  (TF)  at  a  concentration 
of  17  nM  rather  than  the  default  value  of  5  pM.  We  chose 
this  TF  concentration  because,  in  the  average  subject  under 


568  www.anesthesia-analgesia.org 


ANESTHESIA  &  ANALGESIA 


Figure  3.  Model-predicted  quantitative 
parameters  of  thrombin  generation  for 
the  “average”  subject  under  the  normal 
and  hypothermic  temperature  conditions 
with  and  without  dilution.  A,  Clotting 
time  (CT)  (time  to  10  nM  thrombin);  B, 
Thrombin  peak  time  (PT);  C,  Maximum 
slope  of  the  thrombin  curve  (MS);  D, 
Thrombin  peak  height  (PH);  E,  Area  under 
the  thrombin  curve  (AUC);  F,  Prothrombin 
time,  calculated  as  CT  for  initial  tissue 
factor  concentration  =  17  nM.  The  blue 
squares  correspond  to  single  values 
for  the  normal  temperature,  whereas 
the  boxplots  show  thrombin  parameter 
distributions  corresponding  to  groups 
of  thrombin  curves  calculated  for  ran¬ 
domly  generated  sets  of  temperature 
coefficients  (5000  random  sets  were 
generated  independently  for  each  hypo¬ 
thermic  temperature  value  and  each 
dilution  scenario).  In  the  boxplots,  the 
whiskers,  green  horizontal  lines,  and  red 
lines  correspond  to  the  range,  interquar¬ 
tile  range,  and  median  of  the  parameter 
distributions,  respectively.  Statistical 
significance  of  hypothermia-induced  dif¬ 
ferences  (for  each  parameter,  except  PH, 
tested  independently):  normal  tempera¬ 
ture  parameter  values  versus  hypother¬ 
mic  values  for  the  same  parameter,  P  = 
0  to  machine  precision. 
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normal  temperature,  this  TF  concentration  gave  a  CT 
of  11.38  seconds,  which  can  be  considered  a  representa¬ 
tive  value  for  normal  prothrombin  time  in  humans.33  Our 
results  for  thus  estimated  prothrombin  time  (Fig.  3F)  were 
qualitatively  similar  to  those  obtained  for  CT  (Fig.  3A).  Our 
computational  predictions  for  both  prothrombin  time  and 
CT  were  consistent  with  experimental  results  reported  in 
the  literature,7'9,10  including  results  on  combined  effects  of 
hypothermia  and  blood  dilution.8 

To  identify  the  thrombin  generation  parameters  that 
were  affected  by  hypothermia  the  most,  we  calculated  the 
fold  changes  (FCs)  between  the  normal  temperature  values 
of  the  parameters  and  the  corresponding  median  hypo¬ 
thermic  values  (for  the  same  dilution  scenario).  The  FC 
was  defined  for  the  hypothermic  value  (HV)  and  the  nor¬ 
mal  value  (NV)  of  a  parameter  as  follows:  FC  =  FFV/NV  if 
HV  >  NV,  and  FC  =  NV/HV  otherwise.15  For  all  dilution 
scenarios  and  all  temperatures,  the  FC  in  PH  was  1.00.  For 
any  other  parameter  (including  prothrombin  time)  and 
dilution  scenario,  the  FCs  for  31°C,  32°C,  33°C,  34°C,  35°C, 
and  36°C  were  approximately  1.73, 1.58, 1.44, 1.31, 1.20,  and 
1.10,  respectively. 


Temperature  Effects  and  Intersubject  Variability 

Intersubject  differences  in  normal  levels  of  coagulation 
factors  result  in  considerable  intersubject  differences  in 
thrombin  curves  (at  37°C,  Fig.  4A),  which  also  exist  for  hypo¬ 
thermic  temperatures  (e.g.,  at  33°C,  Fig.  4B).  We  thus  inves¬ 
tigated  how  differences  in  coagulation  factor  composition 
impact  the  influence  of  low  temperatures  on  the  quantitative 
parameters  of  thrombin  generation.  To  this  end,  we  applied 
our  simulation  strategy  to  the  group  of  LETS  subjects.  We 
initially  performed  full  TC-randomization  simulations  (with 
randomly  generated  sets  of  5000  TCs)  for  2  LETS  subjects 
having  the  minimum  and  maximum  normal  PT  values  in  the 
group.  These  simulations  demonstrated  that  interquartile 
ranges  for  predicted  thrombin  and  TAT  kinetics  can  be  suffi¬ 
ciently  narrow  even  for  subjects  with  quantitatively  very  dif¬ 
ferent  thrombin  generation  kinetics  (Fig.  4,  C  and  D).  Thus, 
in  our  subsequent  simulations  with  the  LETS  subjects,  we 
focused  on  the  behavior  of  median  thrombin  curves,  which 
we  approximated  by  calculating  thrombin  generation  kinet¬ 
ics  with  each  TC  set  to  2.5  (no  randomization). 

In  each  subject,  decreasing  the  temperature  caused  a 
monotonic  increase  in  CT,  PT,  AUC,  and  prothrombin  time. 
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Figure  4.  Intersubject  variability  in  the 
Leiden  Thrombophilia  Study  (LETS) 
subject  group.  A  and  B,  Thrombin  gen¬ 
eration  in  individual  subjects.  A,  Normal 
temperature  (37°C);  B,  Hypothermic 
temperature  (33°C).  The  green  and  red 
curves  correspond  to  the  subjects  char¬ 
acterized  by  the  smallest  and  the  larg¬ 
est  normal  peak  time,  respectively.  The 
hypothermic  curves  were  calculated 
by  setting  all  temperature  coefficients 
(TCs)  to  2.5  (no  randomization).  C  and 
D,  Temperature  dependence  for  the 
subjects  with  extreme  thrombin  peaks 
in  the  LETS  subject  group.  C,  Thrombin 
curves;  D,  Thrombin-antithrombin  (TAT) 
curves.  Blue  curves  correspond  to  the 
“average”  subject;  green  and  red  curves 
correspond  to  the  LETS  subjects  with 
extreme  thrombin  peaks  (A  and  B).  Solid 
lines  correspond  to  normal  temperature; 
dashed  lines  correspond  to  hypothermic 
temperature  (33°C).  The  upper  and  lower 
dashed  lines  of  the  same  color  desig¬ 
nate  the  interquartile  range  (highlighted 
by  fill  of  the  same  color)  for  a  group  of 
thrombin  curves  calculated  for  5000  ran¬ 
dom  TC  sets. 


and  a  monotonic  decrease  in  MS  (Table  1).  PH  was  practi¬ 
cally  unaffected  by  temperature  in  all  subjects.  Interestingly, 
despite  the  presence  of  noticeable  intersubject  variability  in 
parameter  values  at  each  of  the  temperatures  (Table  1),  the 
FCs  in  the  thrombin  parameters  (and  prothrombin  time)  for 
every  subject  followed  a  quantitative  pattern  very  similar 
to  that  for  the  average  subject  (see  the  Combined  Effects 
of  Temperature  and  Dilution  on  Standard  Quantitative 
Parameters  of  Thrombin  Generation).  Our  results  for  the 
thrombin  curve  parameters  are  consistent  with  analogous 
calculations  for  the  TAT  curve  parameters.  Indeed,  in  our 
LETS  subject  group,  50%  activation  times  for  the  TAT  curves 
were  systematically  increased,  and  maximum  slopes  were 
decreased,  whereas  the  TAT  plateau  was  practically  unaf¬ 
fected  by  hypothermia  (Table  2). 

DISCUSSION 

During  the  initial  phase  of  hemostasis,  blood  from  a  dis¬ 
rupted  vessel  comes  in  contact  with  TF-bearing  cells  in  the 
surrounding  tissues;  this  triggers  the  generation  of  throm¬ 
bin,  which  converts  fibrinogen  into  fibrin,  the  main  struc¬ 
tural  component  of  blood  clots.34-35  Here,  we  extended  the 
validated  Hockin-Mann  kinetic  model17-18  to  systematically 
investigate  how  thrombin  generation  is  affected  by  (sys¬ 
temic  and  possibly  local)  hypothermia,  alone  and  in  com¬ 
bination  with  blood  dilution.  Our  computational  strategy 
was  based  on  the  physically  justified  assumption  that  hypo¬ 
thermia  impacts  thrombin  generation  by  reducing  the  rates 
of  the  biochemical  reactions  in  the  thrombin  generation 
network  (see  Methods  and  Appendix).  Our  results  suggest 
that  thrombin  generation  is  progressively  delayed  at  hypo¬ 
thermic  temperatures  (Fig.  1,  A-D).  While  this  delay  might 
be  slight  (and,  therefore,  hard  to  detect  experimentally)  for 


mild  hypothermia  (34°C-36°C  for  trauma  patients3),  it  is 
noticeably  larger  for  moderate  hypothermia  (32°C-34°C  for 
trauma  patients3)  (Figs.  1-3;  Table  1).  Our  randomization- 
based  computational  strategy  gives  narrowly  localized  pre¬ 
dictions  for  thrombin  generation  kinetics  and  can  allow  one 
to  avoid  TC  randomization  for  certain  analyses  (Fig.  2). 

In  recent  studies,  in  vitro  thrombin  generation  kinetics 
are  often  characterized  by  the  standard  thrombin  curve 
parameters:  CT,  PT,  MS  (the  timing  parameters),  PH,  and 
AUC  (the  amount  parameters).15-16-21  Together,  these  param¬ 
eters  provide  a  more  complete  description  of  all  temporal 
phases  of  thrombin  generation  than  the  traditional  indica¬ 
tors,  such  as  prothrombin  time.  However,  a  comprehensive 
experimental  study  of  hypothermia-induced  effects  on  these 
parameters  has  not  yet  been  reported.  As  a  key  finding  of 
this  study,  our  computational  analysis  of  thrombin  genera¬ 
tion  in  the  average  subject  suggested  that  CT,  PT,  and  AUC 
(as  well  as  prothrombin  time)  are  increased  in  hypothermia, 
whereas  MS  is  decreased,  and  PH  is  practically  unaffected 
(Fig.  3).  Moreover,  we  established  that  the  effects  of  dilution 
on  thrombin  generation  at  hypothermic  temperatures  are 
analogous  to  their  effects  at  normal  temperature  (Fig.  3)  and 
are  fully  consistent  with  the  results  of  our  earlier  study.19 
We  extended  current  knowledge  regarding  multifactorial 
coagulopathy  by  showing  that  hypothermia  and  dilution 
act  in  an  additive  fashion  to  affect  the  parameters  of  throm¬ 
bin  generation  (Fig.  3). 

Our  computations  demonstrated  that,  for  the  average 
subject,  the  magnitude  of  the  hypothermia-induced  FC  in 
every  parameter  affected  by  hypothermia  depends  on  the 
temperature,  but  not  on  the  parameter  itself.  By  applying 
our  computational  strategy  to  the  LETS  subject  group  we 
showed  that,  despite  considerable  intersubject  variability 
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in  absolute  parameter  values  (Fig.  4;  Table  1),  the  relative 
effect  of  hypothermia  on  each  parameter  for  every  subject 
was  nearly  the  same  as  observed  for  the  average  subject. 
While  intersubject  variability  in  hypothermia-induced 
effects  on  thrombin  generation  has  not  been  studied  in  labo¬ 
ratory  experiments,  based  on  our  findings,  we  anticipate 
that  experimentally  measured  relative  effects  of  hypother¬ 
mia  will  only  mildly  depend  on  intersubject  differences.  In 
all  subjects,  the  impact  of  hypothermia  on  thrombin  gen¬ 
eration  predicted  by  our  model  was  substantial.  Indeed,  at 
33°C,  the  CT  and  PT  were  approximately  44%  larger  than 
the  corresponding  normal  values.  This  suggests  that  notice¬ 
able  hypothermia-induced  thrombin  generation  delays  for 
moderate  or  severe  (<32°C  for  trauma  patients3)  hypother¬ 
mia  should  be  expected. 

Our  study  has  limitations  that  follow  from  the  simpli¬ 
fied  representation  of  the  blood  clotting  process  and  our 
choice  of  methodology.  First,  this  is  a  computational  study 
that  does  not  account  for  blood  flow,  as  well  as  for  some  of 
the  components  of  the  blood  coagulation  system,  such  as 
fibrin(ogen),  thrombomodulin,  and  proteins  C  and  S.  Yet, 
this  model  reflects  experimental  setups  that  are  frequently 
used  to  study  thrombin  generation  in  vitro.15'18,19'24  Moreover, 
the  moderate  complexity  of  the  model  allows  for  tractable 
analyses  that  may  elucidate  major  trends  to  be  subsequently 
tested  in  more  complete  (and  complex)  systems.  Second,  our 
modeling  strategy  does  not  explicitly  represent  the  effects  of 
temperature  on  the  cells  (such  as  platelets)  whose  surfaces 
support  thrombin  generation.34  However,  it  is  reasonable  to 
assume  that,  because  many  of  the  underlying  processes  in 
cell  activation  are  governed  by  biochemical  reactions,  the 
dependence  of  cell  activation  on  temperature  can  be  implic¬ 
itly  represented  by  modulating  the  kinetic  constants  in  the 
Hockin-Mann  model  (see  Methods  and  Appendix).  Finally, 
a  lack  of  published  experimental  data  precluded  direct  test¬ 
ing  of  some  of  our  computational  predictions  regarding 
hypothermia-induced  effects  on  thrombin  and  TAT  curve 
parameters.  Particularly,  the  hypothermia-induced  increase 
in  AUC,  predicted  by  our  model  and  experimentally 
detected  for  murine  blood  plasma,36  remains  a  hypothesis  in 
need  of  experimental  validation  for  human  blood.  Being  a 
systematic  study  of  the  effects  of  hypothermia  on  thrombin 
generation,  our  work  sets  a  direction  for  future  experimen¬ 
tal  investigations.  Such  investigations  should  be  aimed  at 
measuring  the  effects  of  in  vitro-induced  hypothermia  on 
thrombin  generation  parameters  in  human  blood  plasma, 
as  well  as  the  effects  of  in  vivo-induced  hypothermia  on 
thrombin  generation  in  porcine  models.  Moreover,  future 
investigations  should  focus  on  the  multifactorial  nature 
of  traumatic  coagulopathy  and  supplement  our  results  on 
hypothermia  and  dilution  with  an  analysis  of  acidosis  and 
acute  traumatic  coagulopathy. 

Our  results  complement  the  existing  experimental  stud¬ 
ies  and  provide  further  evidence  that  thrombin  genera¬ 
tion  impairment  is  a  major  cause  of  hypothermia-induced 
abnormalities  in  blood  coagulation.  A  significant  practical 
consequence  of  this  conclusion  is  the  possibility  to  develop 
pharmacological  strategies  that  reverse  the  effects  of  hypo¬ 
thermia  by  restoring  normal  thrombin  generation.  The 
feasibility  of  this  approach  is  supported  by  the  ability  of 
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Table  2.  Thrombin-Antithrombin  Complex  (TAT)  Curve  Parameters  in  the  Leiden  Thrombophilia  Study  (LETS) 
Subject  Group  Under  Temperature  Variation 

Temperature,  °C 

50%  activation  time,  min 

MS_TAT,  nM/min 

Plateau  level,  nM 

37 

8.27  (7.75-8.88) 

288.76  (249.67-323.14) 

999.71  (901.51-1092.40) 

36 

9.06  (8.49-9.73) 

263.48  (227.81-294.84) 

999.71  (901.51-1092.40) 

35 

9.93  (9.31-10.66) 

240.41  (207.86-269.03) 

999.71  (901.51-1092.40) 

34 

10.88  (10.20-11.68) 

219.36  (189.66-245.47) 

999.71  (901.51-1092.40) 

33 

11.93  (11.18-12.80) 

200.15  (173.06-223.98) 

999.71  (901.51-1092.40) 

32 

13.07  (12.25-14.03) 

182.63  (157.91-204.37) 

999.71  (901.51-1092.40) 

31 

14.32  (13.43-15.38) 

166.64  (144.08-186.48) 

999.71  (901.51-1092.40) 

The  TAT  curves  for  hypothermic  temperatures  (<37°C)  were  computed  with  all  temperature  coefficients  set  to  2.5  (no  randomization).  Fifty  percent  activation  time 
is  the  time  to  50%  of  the  plateau  level.  The  plateau  level  is  estimated  as  the  TAT  level  at  the  end  of  the  considered  time  interval  (50  min).  The  data  are  shown 
as  median  (interquartile  range). 

Statistical  significance  of  hypothermia- induced  differences  (for  50%  activation  time  and  MS_TAT  tested  independently):  normal  temperature  parameter  values 
versus  hypothermic  values  for  the  same  parameter,  P  <  10-5.  The  P-value  computations  were  performed  similarly  to  our  analysis  for  the  data  in  Table  1.  We 
analyzed  1  TAT  curve  parameter  at  a  time;  the  plateau  level  was  insensitive  to  temperature  (see  the  Table)  and  was  therefore  excluded  from  analysis.  For  each  of 
the  remaining  2  TAT  curve  parameters,  we  considered  7  samples,  each  corresponding  to  1  temperature  level  and  comprising  the  472  parameter  values  for  the 
472  LETS  subjects  for  that  temperature.  For  the  7  samples,  we  performed  6  statistical  comparisons  using  Wilcoxon  signed  rank  test:  the  normal  temperature 
sample  was  compared  with  each  of  the  hypothermic  temperature  samples.  Each  such  comparison  returned  a  P  value,  which  were  Bo nferroni -corrected  via 
multiplying  by  6.  Each  of  thus  corrected  P  values  was  <10-5. 

MS_TAT  =  maximum  slope  of  the  TAT  curve. 


recombinant  factor  Vila,  which  acts  at  the  level  of  thrombin 
generation,  to  mitigate  the  adverse  effects  of  hypothermia 
on  blood  coagulation.32'37-38  Yet,  it  is  known  that  other  fac¬ 
tors,  such  as  suggested  impairment  of  platelet  function14'39 
and  reduction  in  fibrinogen  availability,40  may  contribute  to 
hypothermia-induced  coagulopathy.  The  relative  contribu¬ 
tions  and  interdependencies  of  all  these  factors  should  be 
addressed  in  future  investigations.  SB 

APPENDIX 

Supplemental  Methods 

In  mathematical  terms,  the  Hockin-Mann  model  is  a  sys¬ 
tem  of  34  nonlinear  ordinary  differential  equations  that  are 
based  on  the  mass  action  law  of  chemical  kinetics.  We  imple¬ 
mented  the  updated  version17  of  the  model  in  the  SimBiology 
toolbox  of  the  MATLAB  software  suite  (Math Works,  Natick, 
MA).  All  computations  were  performed  in  MATLAB  2012a. 
Our  general  modeling  methodology  followed  our  recent 
work.15'19'20  The  initial  tissue  factor  (TF)  concentration  had 
a  default  standard  value  of  5  pM,  unless  stated  otherwise, 
and  was  not  affected  by  dilution.19-20 

Temperature  coefficients  (TCs)  in  our  modeling  strategy 
were  sampled  uniformly  and  independently  from  the  inter¬ 
val  (2-3),  which  can  be  suggested  as  a  reasonable  choice  of 
sampling  interval.26  TC  randomization  was  deemed  neces¬ 
sary  because  the  simplifying  assumption  of  equal  or  similar 
TCs  throughout  a  biochemical  network  may  not  account  for 
true  system  kinetics.27  Randomization  of  model  parameters 
is  often  used  to  assess  model  behavior  in  the  absence  of 
accurate  estimates  for  kinetic  constants.30-41  Indeed,  deriving 
TCs  from  experimental  measurements  reported  in  the  litera¬ 
ture  was  not  possible  because  the  experimentally  measured 
values  of  the  kinetic  constants  were  adjusted  in  the  process 
of  model  development,18  and  no  experimental  temperature- 
dependence  data  on  these  model-specific  kinetic  constants 
are  currently  available.  The  use  of  a  uniform  sampling  dis¬ 
tribution  is  supported  by  a  statistical  and  information-theo¬ 
retic  principle  (known  as  the  Maximum  Entropy  Principle) 
according  to  which,  in  the  absence  of  detailed  information, 
one  should  use  the  distribution  that  is  characterized  by 
maximum  entropy  under  given  constraints.  In  the  case  of 


a  continuous  distribution  on  a  finite  interval,  the  distribu¬ 
tion  with  maximum  entropy  is  the  uniform  distribution.42 
The  intuition  behind  this  principle  is  that  making  the  most 
"neutral"  choice  of  distribution  may  allow  one  to  minimize 
possible  biases. 

Model  Parameter  Analysis 

To  elucidate  the  model  parameters  that  influence  the  mag¬ 
nitude  of  temperature  effects  the  most,  we  performed 
parameter  sensitivity  analysis  in  the  small  vicinity  of  the 
default  parameter  set  in  the  "average"  subject.  To  this  end, 
we  increased  the  value  of  each  of  the  44  parameters  (i.e., 
the  kinetic  constants)  by  X%  (one  at  a  time),  and  for  each  of 
such  perturbations  we  calculated  |  AFC  | ,  which  is  the  abso¬ 
lute  difference  between  the  fold  change  (FC)  value  returned 
by  the  perturbed  model  and  the  FC  value  for  the  unper¬ 
turbed  model.  The  FC  values  were  computed  for  clotting 
time  (CT),  as  described  in  the  Results  section,  with  hypo¬ 
thermic  temperature  equal  to  32°C.  Here,  we  chose  to  focus 
on  CT  because  hypothermia-induced  change  in  CT  may 
be  regarded  as  a  typical  effect  of  hypothermia,  which  can 
be  detected  both  computationally  and  experimentally.  The 
core  temperature  of  32°C  is  a  critical  temperature  thresh¬ 
old  indicating  the  onset  of  severe  hypothermia  for  trauma 
patients. 

The  resulting  44  values  of  |  AFC  |  (each  correspond¬ 
ing  to  1  perturbed  parameter)  were  sorted  in  descending 
order  to  elucidate  the  most  "sensitive"  parameters.  While 
the  outcome  of  such  ordering  generally  depended  on  the 
value  of  X,  for  X  ranging  from  1  x  10_6%  to  1  x  10~4%,  the 
top  5  most  sensitive  parameters  were  those  with  indices 
9,  4,  10,  2,  and  14  (in  order  of  decreasing  |  AFC  | ).  The  cor¬ 
responding  biochemical  reactions  are  given  in  Table  Al. 
These  reactions  characterize  the  interactions  of  factor  VII 
(FVII)  and  factor  Vila  (FVIIa)  and  their  complexes  with  TF 
and  other  clotting  factors.  These  results  are  consistent  with 
a  previously  performed  sensitivity  analysis  of  the  Hockin- 
Mann  model  at  normal  temperature17  and  indicate  that 
the  reactions  involving  FVII,  FVIIa,  and  TF  may  have  the 
largest  influence  on  the  effects  of  temperature  on  thrombin 
generation. 
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Table  Al.  Biochemical  Reactions  Represented  in 
the  Mathematical  Model  of  Thrombin  Generation17 


Biochemical  reactions 

TF  +  FVII  (1)  <->  (2)  TF:FVII 
TF+FVIIa  (3)  <->  (4)  TF:FVIIa 
TF:FVIIa  +  FVII  ->  (5)  TF:FVIIa  +  FVIIa 
FXa  +  FVII  ->  (6)  FXa  +  FVIIa 
FI  la  +  FVII  ->  (7)  Flla  +  FVIIa 

TF:FVIIa  +  FX  (8)  ^  (9)  TF:FVIIa:FX  -►  (10)  TF:FVIIa:FXa 

TF:FVIIa  +  FXa  (11)  ^  (12)  TF:FVIIa:FXa 

TF:FVIIa  +  FIX  (13)  ^  (14)  TF:FVIIa:FIX  -+  (15)  TF:FVIIa  +  FIXa 

FXa  +  FI  I  -►  (16)  FXa  +  Flla 

Flla  +  FVIII  -►  (17)  Flla  +  FVIIIa 

FVIIIa  +  FIXa  (18)  ^  (19)  FIXa:FVIIIa 

FIXaiFVIIIa  +  FX  (20)  (21)  FIXa:FVIIIa:FX  ->  (22)  FIXaiFVIIIa  +  FXa 

FVIIIa  (23)  ^  (24)  FVIIIal-L  +  FVIIIa2 

FIXa : FVII  la :FX  ->  (25)  FVIIIal-L  +  FVIIIa2  +  FX  +  FIXa 

FIXaiFVIIIa  ->  (25)  FVIIIal-L  +  FVIIIa2  +  FIXa 

Flla  +  FV  — ►  (26)  Flla  +  FVa 

FXa  +  FVa  (27)  ^  (28)  FXaiFVa 

FXa: FVa  +  Fll  (29)  «  (30)  FXaiFVaiFII  -►  (31)  FXaiFVa  +  mFlla 

mFlla  +  FXaiFVa  ->  (32)  Flla  +  FXaiFVa 

FXa  +  TFPI  (33)  <->  (34)  FXaiTFPI 

TFiFVIIaiFXa  +  TFPI  (35)  ^  (36)  TF:FVIIa:FXa:TFPI 

TFiFVIIa  +  FXaiTFPI  -►  (37)  TF: FVI la: FXaiTFPI 

FXa  +  AT  ->  (38)  FXa: AT 

mFlla  +  AT  — >  (39)  mFllaiAT 

FIXa  +  AT  — >  (40)  FIXa: AT 

Flla  +  AT  ->■  (41)  FI  la:  AT 

TFiFVIIa  +  AT  — ►  (42)  TF:FVIIa:AT 

FIXa  +  FX  -►  (43)  FIXa  +  FXa 

mFlla  +  FV  — >  (44)  mFlla  +  FVa 

The  reacting  biochemical  species  are  blood  coagulation  factors  and 
their  inactive  precursors,  whose  description  can  be  found  in  the  original 
publication  that  introduced  the  model18  and  references  cited  therein.  The 
numbers  in  parentheses  pointed  at  by  the  arrows  in  the  reaction  equations 
designate  the  kinetic  constant  indices  for  those  reactions,  e.g.,  A  (x)  <->  (y)  B 
denotes  the  reaction  A  — ►  B  with  rate  y  and  the  reverse  reaction  with  rate  x. 
TF  =  tissue  factor;  AT  =  antithrombin;  TFPi  =  tissue  factor  pathway  inhibitor. 
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